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We review experimental and theoretical tools to excite, study and understand strongly interacting 
Rydberg gases. The focus lies on the excitation of dense ultracold atomic samples close to, or within 
quantum degeneracy, to high lying Rydberg states. The major part is dedicated to highly excited 
S-states of Rubidium, which feature an isotropic van-der-Waals potential. Nevertheless are the setup 
and the methods presented also applicable to other atomic species used in the field of laser cooling 
and atom trapping. 

PACS numbers: 32.80.Ec, 67.85.-d, 42.50.Ct 



INTRODUCTION 

The field of Rydberg atoms has been revolutionized 
by the combination with ultracold atomic gases. The 
availability of commercial high power diode laser systems 
allow for efficient excitation into electronically highly ex- 
cited states with excellent frequency resolution. Here the 
major motivation behind many experiments is to make 
use of the strong interaction mechanism among Rydberg 
atoms, which become especially apparent at large atomic 
densities as found in ultracold gases. On the other hand 
the precise spectroscopy of the interaction energies is 
supported by the availability of trapped atoms, the al- 
most negligible Doppler effect given by small tempera- 
tures and the possibility to prepare all atoms in a spe- 
cific quantum state with high fidelity. Combining this 
control over internal and external degrees of freedom for 
atoms, molecules and ions with switchable long range in- 
teractions establishes an ideal basis for applications in 
quantum simulation, quantum computation and quan- 
tum optics. 

Generally there exist two pathways to prepare a con- 
trollable ensemble of strongly interacting particles. One 
way is to start with a system of weakly interacting par- 
ticles as in the case of quantum degenerate gases and 
then try to implement strong interactions among them 
e.g. with Feshbach resonances [H, optical lattices 
high finesse resonators or by extreme cooling in or- 
der to make the initially weak interaction energy the 
dominant scale Q • The second approach utilizes already 
strongly interacting particles as polar molecules [Hj and 
ions, where the latter already constitutes the best work- 
ing quantum computer Q and universal quantum sim- 
ulator in the world. The drawback of the latter ap- 
proach is, that the presence of interactions requires the 
particles to be restricted to specific geometries. 

At this point Rydberg states might be an ideal add-on 
to both approaches, since they are switchable and exhibit 



interaction strengths comparable to ionic systems 0. 
The major interaction mechanisms between two highly 
excited atoms at large distances is typically given by the 
van-der-Waals interaction and for smaller distances by 
the dipole-dipole interaction . The effects of both in- 
teraction types have been observed in laser cooled clouds 
and also in optically and magnetically trapped 
clouds [H-EH. A prominent effects of the strong inter- 
action mechanism is an excitation blockade in which the 
energy of a Rydberg state is shifted out of resonance with 
respect to the exciting laser by a neighboring Rydberg 
atom. 

The most drastic consequence of the dipole blockade is 
the occurance of a quantum phase transition to a crys- 
talline phase of Rydberg excitations |22l - l24l |. The pa- 
rameters controlling the transition can be tuned with ex- 
ternal laser fields. While this Rydberg crystal has so 
far eluded direct observation, critical fluctuations of this 
quantum phase transition have already been reported (25l | 
in an non adiabatic manner. The realization of the crys- 
talline phase, which is the ground state of the system, 
would require adiabaticity during preparation [26l429j . 

The blockade effect is also used in a more involved 
approach with exactly two atoms located in a certain 
distance establishing the prototype of a quantum gate 
30l 3ll | . Since the loading of individual sites is still prob- 
abilistic [sH, the scalability of this approach is limited 
as long as there is no deterministic loading scheme avail- 
able |33l . |34 1 . An alternative way to achieve a large quan- 
tum register with individually resolved sites is a Mott- 
insulator state with single site resolution [s^ [s^ or a 
switchable local electric field, e.g. produced by an elec- 
tron beam [s^] which allows an addressing via the spatial 
dependent Stark shift. 

Although the first proposal to use the interaction between 
two spatially resolved Rydberg atoms [s^ seemed quite 
involved in the year 2000, it triggered a large theoretical 
effort and several experimental groups started to imple- 
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ment these ideas. There exist elaborate proposals for 

the realization of a universal quantum simulator 

and v arious w ays of to implement a quantum computer 

SIT 



42|-|44|. Some of these schemes rely on collective 



states, where one Rydberg excitation is shared by sev- 



eral atoms [42|, |45|, |46[ . In this article we will focus on 
the many particle character of the collective states in an 
ensemble of ultracold atoms with extensions well larger 
than the blockade radius. A detailed understanding of 
the excitation dynamics in these systems into a strongly 
correlated quantum system, its microscopic description, 
the limitations of the effective models used, the sources 
of dccoherence and dcphasing and finally of the applica- 
ble experimental techniques will not only be helpful for 
the implementation of various quantum simulation and 
quantum computation protocols but also for other tasks 
involving interacting Rydberg atoms. One could be the 
idea of combining interacting Rydberg atoms with quan- 
tum degenerate gases 47H50| to render new interaction 
mechanisms. Besides the effects of the strong interac- 
tion on the material wavefunctions there exists also a 
strong backaction on the radiation field, which can pro- 
duce and also analyze various kinds of non classical light 



fields 26, 51, 5 



In the following we will first describe the basic proper- 
ties of non-interacting Rydberg atoms and deduce from 
this the applicable energy and time scales for meaningful 
experiments. These considerations also set the precondi- 
tions and limitations to the techniques for the excitation 
of ultracold trapped atoms to Rydberg states and their 
subsequent detection. In the next section we will include 
the interaction among Rydberg states and will derive a 
meanfield model for the strongly interacting case, which 
describes quite well the experimental findings. Finally we 
will summarize experimental tools and results on strongly 
interacting quantum systems based on Rydberg atoms. 



tial for the Rydberg electron and a simple 1/r Coloumb 
potential is not applicable anymore. This effect can be 
faced in a phenomenological manner by introducing a 
quantum defect j, /), which results in an altered Ryd- 
berg formula 
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where R' is the specific Rydberg constant for the 
clement in question, e.g. R' = 109736.605 cm^^ for 
Rubidium (R' ^ 109737.316 cm-^ for hydrogen). The 
effective quantum number n* mostly depends on / and is 
in lowest order for Rubidium determined hy S{1 = 0) = 
3.13, S{1 = 1) = 2.64^6(1 = 2) = 1.35, S{1 = 3) = 0.016 
and S{1 > 3) 57H59j. More accurate values for nS 
and nD states have been recently measured [60| . As 
a consequence of the 1-depcndent quantum defects the 
degeneracy of the s,p,d and f-states is lifted and only 
states with I > 3 show a hydrogcnlike behavior. With 
the knowledge of the energies of the individual excited 
states it is pos sible to solve the Schrodingcr equation 
analytically |6l| or numerically to obtain also the radial 
wavefunctions. For that purpose one has to add a term 
Vp = —-^pr to account for the Coulomb potential of 
the core electrons in terms of a core polarizability e.g. 
ad = 9.023 a.u. in the case of Rubdium Since the 
probability of finding a Rydberg electron inside the core 
is quite small the quality of calculated wavefunctions 
are often of sufficient accuracy to be used to reliably 
compute dipole matrix elements, while improved model 
potentials exist for other cases [63|. These dipole matrix 
elements are crucial for the calculation of the excited 
state lifetime, the polarizability, the van-der-Waals 



interaction and many more properties [53 1. 



GENERAL PROPERTIES OF RYDBERG ATOMS 

The properties of Rydberg atoms are similar to that 
of ground state atoms, but are severely enlarged with 
the principal quantum number n (53| . Table H] shows an 
overview of some important values and their correspond- 
ing scaling. 

Electronic structure and lifetimes of highly excited 
Rubidium atoms 

The largest fraction of all laboratories working with 
ultracold atoms use Alkalis, since their simple hydrogen- 
like structure makes their excited states accessible and 
controllable with only a few laser frequencies. For most 
of them the deviation from a hydrogen atom is caused by 
the core electrons which add a repulsive Coulomb poten- 



The lifetime of Rydberg states is determined by the 
radiative decay to lower lying levels but also by transi- 
tions to higher (and lower) lying states induced by black- 
body radiation as in shown in Fig. [T] The major con- 
tribution to the lifetime is given by spontaneous decay 
events to low lying levels where the product of interaction 
strength and density of states of the electromagnetic field 
scaling as lo^ outweighs the decreasing dipole moments 
n*~^-^°..n*~^-^^). For Rydberg state n > AO the ab- 
sorption and stimulated emission of thermally occupied 
infrared modes to neighboring states starts to dominate 
the lifetime. The latter effect can be eliminated by using 
a cryogenic environment as it is done in cavity QED ex- 
periments with Rydberg atoms (65| . 
The transition rate between two states i and / with a 
dipole matrix element {i\r\f), separated by the energy 
huJif, is typically expressed in terms of an Einstein A 
coefficient for the spontaneous decay 
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TABLE L With some explicit values for the 43S state in Rubidium it is possible to get an idea for the framework and the 
boundary conditions for an experimental setting involving Rydberg atoms. With the given scaling laws (n*)^ one can also find 
an estimate for the physical situation with other Rydberg states. 
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FIG. 1: Stacked bar chart for the calculated transition rates 
of the 43S1/2 state of Rubidium (top) and a 30x enlarged 
view (bottom). The Einstein A coefficients for the sponta- 
neous decay are shown in blue nPi/2 and light blue nP-^j^ to 
all energetically lower lying states (5 < n < 42). The sum of 
these coefficients results in a lifetime due to spontaneous de- 
cay of 80.8 \is. An additional loss of population is caused by 
blackbody radiation which favors neighboring states shown 
in red nPi/2 a-nd orange nP-^i^- The total lifetime (includ- 
ing blackbody radiation at 300K) is then reduced to 42.3 /is, 
which corresponds to a Lorentzian linewidth of 3.76 kHz and 
is in good agreement with more elaborate calculations carried 
out in [el. 



(2) 



and an Einstein B coefficient (B = ANiiS)) for stimu- 
lated emission as well as absorption of black body radia- 
tion. The number of black body photons N{uj) per mode 
at a given temperature T is given by 



iV(a;) 
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^huj/kBT _ I ■ 
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which results in a total lifetime t^^ = ^ A + 
J2 AN{lu). The first sum is limited to states with energies 
below the Rydberg state and the second sum runs over 
all allowed dipole transitions, however neglecting transi- 
tions to ionized states. The actual measurement of the 
radiative lifetime by delayed field ionization and detec- 
tion of remaining Rydberg states is quite involved since 
one has to distinguish the Rydberg state in question from 
energetically close by states, which are likely to be pop- 
ulated by blackbody radiation [66[. Here this problem 



has been avoided by a state sensitive detection scheme. 
Another method employs coherent Ramsey spectroscopy, 
which is only sensitive to one specific Rydberg state [67| 
and also delivers reliable lifetimes if the total linewidth 
of the excitation lasers are accordingly small. 



Rubidium Rydberg atoms in electric fields 

The large polarizability of Rydberg states is a welcome 
feature for experimentalists, since it allows for simple and 
fast manipulation of the excited state energies by apply- 
ing small electric fields E. With the knowledge of the 
dipole matrix elements it is straightforward to find the 
new eigenenergies in the presence of an electric field [HI] . 
The corresponding Stark map, as shown in Fig. [2j ex- 
hibits two kinds of Stark shifted states. If the degeneracy 
is lifted by the quantum defect, the Stark effect starts of 
quadratic Estark = o.E'^ jl with the polarizabilty a, and 
is linear in the case of degeneracy, which corresponds 
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1= Ce = (11.97 - 0.8486n + 3.385 x IQ-^n^) . (4) 
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FIG. 2: Calculated Stark map including LS coupling with 
respect to the 43Si/2 state in ^''Rb. Only states with L > 4 
exhibit hydrogen like linear Stark fans. For lower L the states 
are shifted quadratically as it can be seen in the inset for the 
4351/2 state. 



to a permanent dipole moment. In the case of Rubid- 
ium the polarizability for the 43S states is a = —17.7 
MHz/ (V/cm) ^, which shifts this state at a field of only 1 
V/cm already 350 linewidths to the red. For states with 
higher angular momentum the manifold of magnetic sub- 
states discloses a unique pattern in the presence of an 
electric field [g^, which can be used to calibrate the ac- 
tual electric field (zi. 



Various interaction parameters have been calculated 
with pertubative methods fn\ but also by diagonal- 



ization 



Zi-lzil. For Rubidium in the 43S state the 
van der Waals coefficient is —1.697 x 10^^ a.u. which 
is not only much larger than in the ground state 5S 
(Cg = -1-4707 a.u.) but is also opposite in sign resulting 
in a repulsive potential. This is a quite useful feature 
since attractive forces eventually lead to inelastic col- 
lisions and ionization [76l - [80| . At small distances the 
van der Waals interaction becomes comparable to the 
energy differences between the relevant pair states, the 
Forster defect, and pertubative methods arc not valid 
anymore. At smaller distances the interaction potential 
adopts a dipolar ^ character, which is valid down 
to the LeRoy radius where the exchange interaction of 
the electrons cannot be neglected anymore. For the 43S 
state in Rubidium the Forster defect is 3 GHz [HI which 
corresponds to the van-der- Waals interaction energy 
at a distance just below one micron. However, the 
distance between two Rubidium atoms in cvaporativcly 
cooled clouds may easily reach smaller values leading 
to interaction energies of many tens of GHz, an energy 
much larger than all other relevant energies scales in 
the system. Even at a distance of 5 ^m the interaction 
energy is still 150 kHz, which establishes in combination 
with a comparable coupling strength of the driving laser 
fields and sufficiently dense atomic samples the domain 
of strongly interacting Rydberg gases. 



Interaction of Rydberg atoms 

In the case of large distances r between two Rydberg 
atoms the interaction is in most cases determined by 
a van der Waals interaction Vydw = Cq/t^. The van 
der Waals interaction is theoretically best described by a 
second order effect in perturbation theory of the dipole- 
dipole interaction (~ r^^) among two atoms, which leads 
to a scaling ^ r^^ . In zero electric field any atom in 
a single parity eigenstate does not exhibit a permanent 
dipole moment. Nevertheless, fluctuations of the electron 
distribution lead to momentary dipole moments ~ , 
which induce dipole moments in the other atom. In sec- 
ond order perturbation theory the interaction energy is 
proportional to the square of the dipole moment for each 



atom divided by the energy spacing of the pair states [71 



The largest contribution is given by the dipole moments 
to close by Rydberg states, which results in a scaling be- 
havior with respect to the principal quantum number as 
Cg ~ n^^. A numerical evaluation of the interaction po- 
tentials gives e.g. for Rubidium for two Rydberg atoms 
in a S-state and 30 < n < 95 [72| in atomic units a simple 
scaling law 



RYDBERG EXCITATION OF ULTRACOLD 
ATOMIC GASES 

For experiments with Rydberg atoms in ultracold gases 
additional equipment and methods for excitation and de- 
tection of Rydberg states are required, without affecting 
the cooling and trapping procedure. The required level 
of performance of these new tools is determined by the 
properties of the Rydberg atoms. For example the laser 
light for the Rydberg excitation has to be sufficiently 
intense and narrow in frequency to achieve a coherent 
evolution within the natural lifetime. The 43S in Rubid- 
ium has a lifetime (including BBR at 300K) of 42.3 ^s 
which corresponds to a linewidth of 3.76 kHz and the ex- 
citation laser should be accordingly narrow. The lifetime 
also sets the timescale for the experiments and by this the 
desired coupling strength of the laser system. To observe 
e.g. a full Rabi-oscillation within the radiative lifetime 
of the 43S state in Rubidium, the Rabi frequency VLr has 
to be accordingly strong and the linewidth of the exci- 
tation sufficiently narrow. Given that the dipole matrix 
elements for excitation into Rydberg states are typically 
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quite small, a sufficient strong laser source with several 
tens of milliwatts is appropriate. But many proposals 
e.g. for quantum computing and quantum simulation 
will only be manageable with several hundred milliwatts 
of available laser power. 

After excitation one also wants to detect the produced 
Rydberg atoms efficiently. Absorptive and fluorescing 
methods, as typically used for ground state atoms, are 
impractical since the spontaneous scattering rates are 
roughly a thousand times smaller. One exception are 
Rydberg states of alkaline earths, where the second non- 
excited electron can be used for optical detection [s^ . 
The most common method relies on efficient and fast 
ion detection with mutlichannel plates or channeltrons, 
which can reach counting efhciencies of 90%. To do so 
one has to field ionize the Rydberg states first cither with 
a strong electric field or by a far infrared laser pulse [sl] . 
The field plates inside the vacuum chamber to ionize the 
atoms and to guide them to the ion detectors are also 
very useful to apply moderate electric fields to shift the 
Rydberg levels in and out of resonance with the driving 
laser fields, to subtract unwanted ions or to implement a 
Ramsey interferometer [84 1. 

The other contributor, ultracold atomic gases, arc 15 
years after the fist observation of Bose-Einstein conden- 
sation [s^ [s^ of the shelf in many atomic physics labo- 
ratories. There exist countless tools and methods to con- 
trol the internal and external degrees of freedom, where 
most of them also get along with the technical constraints 
of highly excited Rydberg states. For precision spec- 
troscopy trapped gases at extremely low temperatures 
arc an ideal sample. The Doppler-effect is largely re- 
duced, the interaction time with the probing light field 
is limited either by the excited state lifetime or by the 
lifetime of the trap, which can be minutes. Actually on 
the typical timescales of Rydberg experiments of a few 
microseconds the ultracold atoms move only a fraction of 
the wavelength of the exciting lasers, which has led to the 
expression of frozen Rydberg gases. Additionally to this 
the high densities of fully polarized samples make the in- 
teraction potential among Rydberg states accessible to a 
spectroscopic precision limited only by the excited state 
lifetime. Such measurements can also be taken in zero 
electric and magnetic fields by switching of all trapping 
potentials (magnetic or dipole) during excitation. 



Two-photon excitation of ultracold atoms into 
Rydberg states 

The excitation into Rydberg states with principal 
quantum numbers ranging from rt = 20 up to the ioniza- 
tion threshold is usually accomplished by a two-photon- 
and sometimes by three-photon excitation schemes. For 
a single-photon excitation of Rubidium to Rydberg states 
one would have to use light at 297 nm, which is difficult 
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FIG. 3: The left part of the picture depicts the level scheme of 
Rubidium, which is well described by the quantum defects as 
described in the text. Due to angular momentum conservation 
it is only possible to reach s- and d-states with two photons 
when starting from the 5s ground state. In all experiments we 
excite the atoms from the 5s F=2,mF=2 via the 5p F=3 mf=3 
to the 43s j=l/2 state. Although the hyperfine splitting of the 
43s state is roughly 500 kHz [H^, HOl we do not observed the 
F' — 1 line due to the selection rules of the given excitation 
scheme. This three-level system can now be reduced to an 
effective two level system by introducing a large detuning S, 
which largely simplifies the subsequent modeling. 



to produce and would also not allow to access the de- 
sired S-statcs. In Rubidium there are in principle two 
two-photon excitation schemes feasible, here we use light 
at 780 nm and 480 nm for to the 5s-5p-ns/nd scheme 
as shown in Fig. [3] but also the alternative route 5s- 
6p-ns/nd with 420 nm and 1016 nm lies in the range of 
standard diode laser systems. 

In all experiments, we detune the red laser far from 
resonance to minimize spontaneous scattering of photons 
via population of the 5p state, which has a lifetime of 26 
ns [87| . In doing so we choose a blue detuning to avoid all 
the lower lying hyperfine levels of the 5p state as shown 
in Fig. [3l For large enough detunings S the intermediate 
state can be eliminated and one ends up with an effective 
two level system with a total coupling strength rin = 
y^TjpT^/ifP + A^ . To obtain coherent coupling between 
the ground state and the Rydberg state it is necessary 
to achieve an effective Rabi-frequency f2 well larger than 
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the linewidth of the Rydberg state (3.76 kHz for 43 S) 
or of the driving laser fields. This is easily achieved for 
the 5s-5p transition with only a few milliwatts of laser 'n' 
power but the dipole matrix element for the second step 5 
(5p-43s) is roughly 10^ times smaller requiring 10^ more £ 
laser power for a comparable coupling strength. 

The typical excitation scheme excites in a first step a 
spin polarized sample in the 5S1/2F = 2,mF = 2 state 
into the 5P3/2F = 3,mp = 3 state and in the second 
step into the n5'i/2 J = 1/2, = 1/2 Rydberg state. 
Note that the hyperfine structure of highly excited state 
can not be resolved and we return here to the J, nij basis. 
Following the standard notation 88, 8^ we can compute 
the dipole matrix element for the second step by 



{Fmf\erq\J'm,p} = {J\\er\\J') x 



(5) 



J / F \ I' J' 1 J 

mji + q mp — mji — q —mp J \mji q —mji~ 

The reduced matrix clement ( J|| er|| J') is related to the 
pure radial part {RnL\er\Rn'L') ~ j RnLfRn'L'T^dr as 

{J\\er\\J') = {RnL\er\R,,.L') x (6) 
V(2J' + 1){2L + i)(_i)./'+^+s+i+(i'-L+i)/2 X 



2i + l 



L L' 1 
J' J S 



For the transition form the 5P3/2 into the 43Si/2 
state we obtain the purely radial matrix element 
(i?5_i|er|i?43^o) = 0.01786 eag leading to a total dipole 
matrix element of 
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FIG. 4: The Doppler broadening depends on the directions of 
the incident laser beams. For parallel irradiation the Doppler 
shift of the two wavelengths adds up, whereas in the counter 
propagating case it is partially compensated. The Doppler 
width has to be added to the lifetime limited Lorenzian hne- 
shape of typically narrow Rydberg lines (e.g. 3.8 kHz for the 
43S state) in terms of Voigt profiles. 



the Gaussian widths Av ~ ■\/8fctr In 2/ mc^ {vi, ± Vr) of 
a Doppler broadened spectrum for a colinear (=^) and a 
counter propagating (<=^) excitation are shown, which has 
to be convoluted with the T = QK spectrum. Here hi/r 
and hiyh are the energy spacing of the lower transition 
and respectively the upper transition. 



Experimental setup for frozen Rydberg gases 

Laser system 



{5P3/2,F^2,mp 
= O.OlOSleao 



2|er,=+i|435i/2, J' = 1/2, mj, = 1/2) 

(7) 



State of the art diode laser systems based on tapered 
amplifiers with subsequent frequency doubling deliver 
about P = 300 mW of blue light, which corresponds to an 
intensity of / = 9.55 MW/m^ at the center of a focused 
Gaussian beam with a width of wq = lOOfim. With the 
electric field E = ^21 je^c we can compute the coupling 
strength i7b ~ dE/h — 11.2 MHz. In combination with 
a coupling strength for the first transition of Vtpi = 100 
MHz and a detuning 6 = 400MHz we obtain a resonant 
(A = 0) two photon coupling strength = 1.4MHz. 

For higher intensities one can consider dye jet lasers 
or dye amplifiers based on Bethune cells where the lat- 
ter has been shown to produce coupling strengths in the 
GHz regime 

Another important consequence of rather narrow 
linewidths of Rydberg states is that the Doppler broad- 
ening becomes comparable to the temperature of ultra- 
cold atomic clouds in the microkelvin regime. In Fig. 2] 



As already mentioned above we excite Rydberg states 
with principal quantum numbers ranging from n = 20 up 
to the ionization threshold by a two-photon excitation 
scheme. The lower transition from the 5S ground state 
to an intermediate 5P state is accomplished with a 
wavelength at 780 nm and the second step to the desired 
Rydberg state with roughly 480 nm. Light at 780 nm 
and 960 nm is produced by standard extended cavity 
diode laser setups in master slave configuration. The 
slave lasers are seeded via acousto-optical modulators 
(AOM) in double-pass configuration which allows to 
tune the light frequency over a range of 60 MHz and 
320 MHz, respectively, within a few microseconds. 
The infrared light is further amplified using a tapered 
amplifier and then frequency doubled (TA-SHG 110, 
Toptica Photonics AG, Germany) to generate the needed 
light at 480nm. 

Note that the experimental results on strongly in- 
teracting Rydberg gases Igil - fgsl l have been 
accomplished with a precursor to our current setup 
based on a transfer cavity with a total linewidth of <1 
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MHz. With the new setup shown in Fig. [5] we were 
able to reduce the hnewidth of the combined two step 
excitation scheme to 60 kHz, which is an improvement 
of more than one order of magnitude compared to 
the old setup. The key element of the new setup are 
two passively stable reference cavities 9J] on which 
we lock both lasers respectively. Therefore we use a 
Pound-Drever Hall-scheme [9^1 with a fast feedback 
(bandwidth«10 MHz) on the current of the laser diode 
and a slow feedback on the grating of the extended 
cavity. To maintain a preferably constant length of the 
plano-concave cavities, the mirrors with a reflectivity 
of 0.995 were mounted on one single Ccrodur spacer 
with a length of 100 mm which corresponds to a free 
spectral range of 1.50 GHz and a Finesse of 626. To 
avoid thermal drifts we mounted the cavity inside a 
high vacuum chamber, vibrationally isolated by nitrile 
o-rings, at a pressure better than 10 mbar. Finally the 
whole setup is set into an actively temperature stabilized 
case. One concave mirror is additionally equipped with 
piezos to tune the length of the resonator. This setup is 
suitable to hit all the desired wavelengths for arbitrary 
Rydberg states, but may be bypassed by the usage of a 
more broad-band AOM-based setup, the sidebands of 
widely tunable EOMs, or by locking to higher TE modes 
of the cavities. The latter comes along with the problem 
that each TE mode exhibits a slightly different lineshape 
requiring adapted parameters for the FID locking circuit. 
By the special arrangement of the piezos, which exhibit 
a 100 times larger temperature coefficient than Cerodur, 
shown in Fig. [6] one can compensate reasonably well for 
thermal drifts. An even more advanced setup would 
involve ultra-low expansion glasses (ULE) as a spacer 
to reduce the thermal drifts even further. Additionally 
we use very stable high voltage supplies (TIDP 005 
106 EPU; iseg Spezialelektronik GmbH, Germany) with 
a voltage stability below 10~^. Since the piezos only 
approach their equilibrium length exponentially after a 
sudden change of the applied voltage we have to wait 
for a few hours to be within the stability range of our 
locking scheme. The combination of all these elements 
result in a frequency stability of lOOkHz/mK, which can 
be further improved by the usage of ultra-low expansion 
glass (ULE) and leaving out the Piezo-rings. 
During the experiment we can tune the frequency by the 
double pass AOMs between the master and slave lasers 
which results in an overall scan range of 640 MHz of the 
blue and 60 MHz of red light. 

The short time stability of the laser system has been 
studied by measuring the linewidth A/ of the light emit- 
ted by the red and infrared slave lasers respectively us- 
ing a delayed self- heterodyne interferometer (DSHI) [o^ . 
The setup is depicted in the inset of Fig. [71 In one path, 
the laser light is shifted in frequency using an AOM at 80 
MHz; in the other arm the light is delayed using a single- 
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FIG. 5: Laser-system for two photon excitation of *^Rb into 
Rydberg states. The red light at 780 nm is generated by a 
standard external cavity diode laser system. The production 
of blue light in the range of 475 nm to 483 nm is more in- 
volved. For this purpose we use a master slave setup, where 
a standard diode laser setup delivers infrared light at 960 nm 
which is amplified by a tapered amplifier. A subsequent fre- 
quency doubling cavity delivers the desired wavelength. Both 
lasers are stabilized (780 nm and 960 nm) to a passively stable 
resonator. 
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FIG. 6: The stability of the laser system is derived from the 
passive length stability of a Zerodur-glass spacer of a Fabry- 
Perot interferometer (FPI). To avoid thermal drifts the whole 
FPI is placed inside a UHV-chamber which is again place 
in a temperature stabilized environment. The length of the 
plan-concave resonator is 100mm and the finesse has been 
measured to F = 626. The spacer contains in total four FPIs 
where two of them are additionally equipped with Piezo rings 
for fine tuning of the length. To minimize thermal drifts the 
Piezo rings are mounted in a self compensating way. 



mode optical fiber of 10.5 km corresponding to a delay 
time of r = 52ps . For coherence times 1/A/ shorter 
than the delay time r, the two beams are completely 
decorrelated. The delayed light can then be treated as 
a second independent but otherwise identical laser. The 
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FIG. 7: Spectral power density of the beat signal from the 
slave laser at 780 nm with master laser locked (blue) and 
free running (red), measured with a delayed self-heterodyne 
interferometer (DSHI) as depicted in the inset. The curves 
are centered around the frequency of the AOM at 80 MHz. 
The signal of the free running laser shows a Voigt-profile with 
a fitted linewidth of 210 kHz (black line) whereas the curve of 
the locked laser shows typical features occurring if the laser 
linewidth is well below the resolution limit of the DSHI (here 
« 19 kHz). 



spectrum Sj of the beat signal from the DSHI is then 
given as the cross correlation of the original laser spec- 
trum Se |97|: 



Si{uj + uj')Si{uj')duj' 



(8) 



By assuming a model for the laser spectrum Si this 
allows us to determine the linewidth of the laser quanti- 
tatively down to the resolution limit 1/t = 19 kHz, given 
by the delay time. The timescale of this measurement 
is also determined by the delay time r = 52p.s which 
is adapted to the typical timescale of our experiments. 
For linewidths below the resolution limit, the spectrum 
of the beat signal shows typical features such as a well 
pronounced peak at the center and a modulation at 
the wings whose periodicity depends on the delay time 
T only. In this range only qualitatively measurements 
are possible unless the laser shows a pure Lorentzian 
spectrum [98| . Both diode lasers in our new setup show 
such features. Considering the frequency doubling and 
the residual broadening of AOMs used for switching 
we can state that combined linewidth of our Rydberg 
excitation laser system is below 60 kHz, which agrees 
with the observed linewidths of experimental Rydberg 
spectra. 

In the end the light is brought to the atoms by po- 
larization maintaining fibers and overlapped by dichroic 
mirrors. At the position of the atoms the 1/e^ radius of 



the 780 nm light was set to 550 pm and the one at 480 
nm to 35 pm. The maximum available laser power in the 
blue light (480 nm) is about 100 mW at the position of 
the atoms. In most experiments the power of the 780 nm 
light is reduced well below one milliwatt to avoid exci- 
tation into the 5P3/2 state. At a detuning of 400 MHz 
(see Fig. [11]) and a typical laser power of e.g. 50 pW the 
spontaneous scattering rate reduces to below one 1 kHz. 
The effective two photon Rabi frequency at this setting 
is 250 kHz. Another important aspect is the uniformity 
of the illumination of the atoms. An atomic cloud at e.g. 
3.4 pK confined in our magnetic trap at an offset field 
of 0.89 G, as used in [l^, has a Gaussian shape with a 
radial width of Up = 8.6 pm. At this parameters 85 % of 
the atoms experience at least 80 % of the maximum two 
photon Rabi frequency. 



Main vacuum chamber 

Our vacuum chamber for exciting and detecting 
Rydberg atoms in an ultracold cloud of atoms is inspired 
by the design of the MIT group described in The 
vacuum chamber consists mainly of two parts, where the 
main chamber is shown in Fig. |8l which includes now 
field plates to apply various electric fields and multi- 
channel plates for detection of charged particles. Not 
shown is the effusive oven assembly, which is operated 
at high vacuum (10~^ mbar) and delivers a thermal 
beam of gaseous Rubidium atoms to an increasing field 
Zeeman slower. The precooled atoms are loaded into a 
magneto-optical trap in the ultra-high vacuum part of 
the setup (< 2 • 10^^^ mbar). After a short molasses we 
transfer the atoms into a magnetic trapping potential. 
By forced evaporative cooling we produce for most 
experiments described here ultracold samples of several 
10^ atoms at a few /xK in the F = 2, mp = 2 state. This 
is done to obtain a pure Gaussian density distribution 
in the harmonic trapping potential, which simplifies the 
subsequent analysis of our data. Typically, we work 
at peak densities of up to 10^^ cm"'^ which is about 
four orders of magnitude larger than in experiments 
with laser cooled samples. If desired we can also extend 
our evaporative cooling to quantum degeneracy and 
produce Bose-Einstein condensates with a few 10^ 
atoms. Note that for the excitation of Rydberg states 
in a Bose-Einstein condensate one has to include in 
addition the bimodal density distribution of partially 
condensed clouds into the analysis as it has been done 
in liil. 

The main chamber depicted in Fig. [8] has to fulfill 
several boundary conditions simultaneously. First of all 
a good optical access in three dimensions for laser cooling 
is mandatory. Two further optical axes are added for 
imaging which are equipped with larger view-ports (CF 
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63) to obtain a better optical resolution. The numerical 
apertures in both imaging axes are 0.17 which yields an 
optimal resolution of 5.6 pm for a wavelength of 780 
nm. Altogether eleven optical view-ports are available 
which are suitable for wavelengths ranging from 300 
nm to 2.5 pm. At the same time close by magnetic 
coils for magnetic trapping are necessa ry. We use a 
cloverleaf style lofFe-Pritchard trap [9^ Il00l |. which is 
located inside two recessed bucket windows outside the 
vacuum. The inner spacing of the two coil assemblies 
is 32 mm. The two pinch coils of the cloverleaf trap 
produce an axial curvature of B" = 0.56 G/cm^ per 
Ampere and the leaves a radial gradient of B' — 0.61 
G/cm per Ampere. At typical operation conditions 
with 400 A in all coils and an offset field of 1.5 G wc 
obtain for ^^Rb trapped in the F = 2, mi? = 2 state 
trapping frequencies of 18 Hz axially and 250 Hz radially. 

The main goal of the setup described here is the inves- 
tigation of Rydberg atoms excited from ultracold gases 
or Bose-Einstein condensates. For the manipulation and 
detection of the highly excited atoms wc included eight 
field plates and two multi-channel plates (MCP) close 
to the atoms inside the vacuum. The details of this 
add-ons will be discussed in detail below. Because of 
the arrangement of the field plates we had to relocate 
the radio-frequency coils used for evaporative cooling 
further away from the atoms as can be seen in Fig. [5] 
and Fig. [S] We use two coils made of polyimide coated 
copper wires consisting each of two loops. Remanent 
charge on the insulating coating can cause disturbing 
electric fields if they arc too close to the atoms, which 
is also avoided by the larger distance of the coils. 
Nevertheless, the coils produce at the position of the 
atoms an average magnetic field of 5 mG (for frequencies 
from 1 to 30 MHz), when driven with 2 W at 50 Ohm 
without impedance matching. This is sufficient to drive 
the magnetic dipole transitions for evaporative cooling. 

A very useful addition to the setup is a helical an- 
tenna operated from outside the vacuum chamber and 
optimized for 6.8 GHz, which corresponds to the hyper- 
fine splitting of ^^Rb. With this antenna it is possible 
to adjust the density of the F = 2, mp = 2 atoms by 
transferring a certain fraction by a Landau-Zener sweep 
to the magnetically untrapped F — \ mp = \ state. Dur- 
ing this process we do not observe a severe change of the 
temperature, which implies also an invariant density dis- 
tribution. By this we avoid the influence of all size effects 
during Rydberg excitation. 



Electric field manipulation and detection of Rydberg atoms 

The high sensitivity of Rydberg atoms to electric fields 
opens the possibility to manipulate the internal states 



of the Rydberg atoms by field plates 10ll | . To produce 
electric field configurations as versatile as possible, we 
installed eight field plates close to the atoms. The 
spatial arrangement can be seen in Fig. |8] and Fig. HI 
Each of these plates can be addressed individually, which 
allows us to generate various field configuration. For 
example for a constant electric field of 1 V/cm pointing 
along the z-axis one has to charge the plates A — D 
to -2.6 V and the plates E — H to +2.6 V/cm, which 
produces in addition a curvature below 0.5 V/cm^ in all 
directions. This corresponds to an additional electric 
field of 50 mV/cm at a displacement of 1 mm with 
respect to the geometric center. Actually we observe 
such an offset field around this order of magnitude [20| . 
which can be explained by a spatial mismatch of the 
magnetic trapping potential with respect to the electric 
field plates. In principle it is possible to compensate 
for such additional fields, but it requires an elaborate 
calibration procedure in all three dimensions, since the 
direction of the perturbing field is most likely unknown. 

The eight field plates are made of stainless steel with 
a thickness of 0.5 mm. They are glued (Epotek 377, 
Epoxy Technology) with 1 mm thick ceramic spacers 
onto the recessed buckets. The dimensions of the spacers 
have to be small enough, that they are completely 
hidden behind the field plates from the viewpoint of the 
atoms. Any insulating surface can accumulate charge 
and falsify the desired field configuration. To charge the 
field plates they are spot welded to stainless steel wires, 
which are radially led outwards as can be seen in Fig. El 
At the edge of the recessed bucket the wires are fixed in 
position each by a short ceramic tubing, which is also 
glued to the buckets and are subsequently connected 
to capton insulated copper wires. These copper wires 
arc then finally connected to one of the fourfold high 
voltage feedthroughs. To avert breakthroughs inside the 
chamber induced by sharp edges we rounded off all four 
edges of each plate with a radius of 1.5 mm. Finally, 
we etched and electro-polished all field plates including 
the spot welded wires to burnish also small spikes. The 
polishing was done in a acid bath consisting of one part 
of 96% sulfuric acid, two parts of 85% phosphoric acid 
and six parts of distilled water. After two minutes at 
a current of 5 A about 70 pm of stainless steel from 
the plates was removed and they exhibited a semi 
gloss surface. After installation of the field plates and 
evacuating the chamber we measured no current leakage 
up to 3000 Volts for all plates. 

During the experimental process it is necessary to 
switch the applied voltages within short times. To do so 
we use bipolar high voltage switches (HTS-6103 GSM, 
Behlke Electronic GmbH, Germany), which have an in- 
trinsic rise-time of 60 ns. The push-pull circuit of the 
switch has to be adjusted to match the capacitive load 
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FIG. 8: Detailed view of the main chamber. In the left part of the figure one recessed bucket is removed to reveal the 
components used for electrical field manipulation and Rydberg atom detection. The right part of the figure shows a slice of the 
main chamber including both recessed buckets. Inside the buckets, but outside the vacuum, the coils for magnetic trapping are 
located. More details are given in the text. 




FIG. 9: Electric field plates (A-H) and Faraday cages (/ and 
J) for the multi-channel plates. Four field plates are glued 
onto each of the recessed buckets, such that plate A, B, C, 
and D lie vis-a-vis to the plates E, F, G and H. The inner 
distance between the plates is 25 mm. All dimensions given 
in the figure are in millimeters. The MCPs were located as 
close as possible to the center of the vacuum chamber without 
loosing any optical access, which resulted in the two different 
distances. 



of 50 pF of each field plate as well as the 300 pF load of 
the high voltage coax cable, which connects the switch 
to the high voltage feedthroughs. 

For a high detection sensitivity of Rydberg atoms, we 
installed two MCPs (Type B012VA, El-Mul Technologies 
Ltd, Israel) inside the vacuum chamber. After field ion- 
ization of the Rydberg atoms, we use one MCP to detect 
the ions. The second MCP is designated to detect simul- 



taneously the electrons. To improve the amplification 
even further we use MCPs in a Chevron configuration, 
which consist of two successive glass plates with a small 
spacing in-between. The electron current arriving at the 
anode is converted by a large resistor to a voltage and 
then amplified by a homebuild circuit including a low 
noise operation amplifier. The whole MCP Chevron as- 
sembly is boxed into a Faraday cage in order to shield 
the atoms in the center of the chamber from the biased 
front side, typically charged with -2000 V. The Faraday 
cage is closed at the front by a grid with a diameter of 
12 mm and a transmittance of 85%. The active area of 
the MCP front side has a diameter of 8.5 mm. 

To detect the Rydberg atoms with a MCP one has first 
to field ionize the excited atoms. This can be done by a 
large enough electric field, which is in the case of a 43Si/2 
state about 160 V/cm. In our case we want to detect 
the ions and one has to provide besides a sufficient field 
strength also a suitable electric field distribution which 
guides the ions into the upper MCP. 

Usually the magnetic fields of the trapping potential 
are still switched on when the ions move towards the 
MCP. The combination of electric and magnetic fields 
provoke a drift on the ions according to the force F = 
q{E + -i/ X B). An estimation for the given experimental 
situation shows that the drift is in our case of the order 
of 1 mm and by this well below the aperture of the MCP 
of 8.5 mm. 

We calibrated the MCP ion signal by monitoring the 
losses in a cold atomic cloud due to Rydberg excitation 
and the corresponding voltage signal on the anode. Af- 
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FIG. 10: Typical experimental sequence. In a first step we 
prepare within 40s a magnetically trapped cloud of ultracold 
atoms. During the two photon excitation we apply a small 
electric field to extract unwanted ions from the sample. After 
excitation we raise the electric field to 200 V/cm and detect 
the field ionized atoms with an MCP detector. The excitation 
and detection of Rydberg states can be repeated up to 600 
times within one atomic cloud. Finally we take an absorption 
image of the remaining atoms. 



tcr amplification, which is the same for the subsequent 
experiments, we acquire a signal of 1 Vs per 3.65 • 10^*^ 
atoms. In principle one could distinguish between sin- 
gle ion events, but the noise level of the signal limits our 
minimum sensitivity to about one hundred ions. Since 
we shoot around 10^2 ions per year into the detector we 
do observe aging effects of the MCP, which makes a re- 
calibration from time to time necessary. Due to the large 
distance of the MCP to the position of the atom cloud 
we can easily distinguish between Rb"*" and RbJ ions by 
their different time of flights due to their differing masses 

fini. 



Excitation of Rydberg states in a magnetic trap 

Due to the confinement of the cold atomic cloud in 
a magnetic trap Rydberg atoms are excited in an off- 
set magnetic field. The experimental sequence for Ryd- 
berg excitation in a magnetic trap is shown in Fig. 1101 
The ultracold cloud consists typically of 4 • 10^ atoms 
at a temperature of 3 fiK confined in a cigar shaped 
harmonic trapping potential with an axial trapping fre- 
quency (along the z-axis) of 18 Hz and an radial trapping 
frequency of 310 Hz. The magnetic offset field at the cen- 
ter of the trapping potential was set to 1.0 Gauss during 
the cooling and trapping sequence. The cloud is spin po- 
larized in the F = 2,mp = 2 ground state with respect 
to the quantization axis given by the magnetic field af- 
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FIG. 11: Level scheme of the relevant levels 5Si/2, 5Pi/2 and 
43Si/2- The left side shows a schematic of the two-photon 
transition including the frequency splittings (in MHz unless 
otherwise noted) between the hyperfine states (nuclear spin 
1=3/2). The right side shows the magnetic sublevels together 
with their respective Zeeman shift in f-basis for 5Si/2, 5Pi/2 
and in j-basis for 43Si/2. The multiple lines denote the num- 
ber of states of the respective other basis each state is com- 
posed of. 



ter reloading from the MOT to the magnetic trap. Since 
evaporative cooling is performed via rf-coupling of neigh- 
boring mp states, a residual occupation of the weakly 
bound F = 2,mp = 1 state is created in the cooling se- 
quence. 

For highly excited states the hyperfine coupling is signif- 
icantly decreased and we remain with the coupling of the 
electron spin and the orbital momentum to a total mo- 
mentum J as shown in figure [TT] Atoms are excited to 
the Rydberg state via a two-photon excitation, detuned 
by about 400 MHz from the intermediate 5Pi/2 state (red 
and blue solid arrow in figure [TT]). The nS 1/2 Rydberg 
states exhibit the same magnetic moment as the ground 
state. Therefore the transition from the ground state to 
this state is insensitive to magnetic fields. Polarizations 
of both excitation lasers are chosen as indicated in figure 
llll to suppress excitation via different paths (gray dashed 
lines) to other magnetic substates. However, the direc- 
tion of the magnetic field is spatially inhomogeneous due 
to the geometry of the magnetic trap. This causes spa- 
tially dependent excitations into the other mj Rydberg 
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states by admixtures of other polarizations. 
Atoms in the F = 2, nip = 1 ground state can also be 
excited to the same Rydberg states. Due to the Zeeman 
splitting of the ground state these transitions show differ- 
ent dependencies on the magnetic field. To determine the 
relevant transitions and explore the dependence on the 
magnetic offset field, Rydberg spectra at different mag- 
netic field strength are taken as shown in Fig. [T2|) . The 
magnetic offset field was ramped up shortly before the ex- 
citation and detection sequence. The spectra show four 
distinct resonances of which three are dependent on the 
magnetic field (solid lines in Fig. [T2|) . These lines can be 
assigned to the expected transitions from F ^ 2, mp = 2 
and F = 2,mF = 1 to 43Si/2,'mj = ±1/2 by their 
magnetic field dependence. The F = 2,mF ~ 2 to 
43S'i/2,'7ij = -1-1/2 transition, labeled by (3) in Fig. [T^ 
is the desired excitation path. 



STRONGLY INTERACTING RYDBERG GASES 

The strength of interactions in a system is usually char- 
acterized by the dimensionless quantity na^ , with n being 
the density and a an effective range of the interaction po- 
tential. Strongly interacting systems with na'^ 3> 1 are 
often difficult to describe theoretically. This can be un- 
derstood as strong interactions typically lead to strong 
quantum correlations, where an effective description in 
terms of non-interacting quasiparticle excitations is no 
longer valid. The first steps to understanding the es- 
sential properties of a strongly interacting systems are 
usually based on mean-field theory, however, especially 
in low-dimensional systems this can lead to unreliable 
results. Nevertheless, it is possible to derive a mean-field 
approach to strongly interacting Rydberg gases that re- 
produces the correct critical exponents for the underly- 
ing quantum phase transition if the dipole blockade is 
explicitely taken into account 22, l25| . In the case of in- 
teracting Rydberg gases the effective range of the van dcr 
Waals interaction between Rydberg states is on the order 
of a « 5/Ltm, which has to be compared with the inter- 
particle distance between the atoms in a Bose-Einstein 
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which is deep in the strongly interacting regime. 
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Interacting Rydberg Gases 
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FIG. 12: Excitation of Rydberg atoms into the 435*1/2 state 
from a cloud of magnetically trapped atoms with a tempera- 
ture of 3 i-iK. a) shows the transitions from the two trapped 
magnetic substates of the ground state to the two magnetic 
substates of the Rydberg level. Transition 3 is the desired ex- 
citation path that is magnetic field independent, b) shows 
spectra of Rydberg excitation at different magnetic offset 
fields ranging equidistant from 0.86 G to 9.34 G in a wa- 
terfall plot. Four resonances are visible whereof three are 
dependent on the magnetic field (see text). Solid lines show 
a linear fit to the peak centers, consistent with the expected 
Zeeman-shifts. Note that the ion detector is saturated for 
the strong lines in the spectra. Therefore a comparison of the 
transition strength is not possible. The resonance at -5.5 MHz 
marked by the dashed line is due to ultralong-range Rydberg 
molecules llOSH. 



The interaction energy per particle in a gas of highly 
excited Rydberg atoms with density n is given in 
the case of a van-der- Waals- interaction by Cqu^. For 
rubidium atoms in the 435' state Cg reaches — 1.7 x 10^^ 
a.u. = -hx 2441 MHz /im®, which corresponds with 
a typical density of n = 10^^ cm~'^ (mean distance 
l/v^ = 1/im) in an ultracold atomic cloud to a pairwisc 



interaction energy of —h x 2441 MHz. This energy 
has to be compared to the other energy scales of the 
system like the temperature in the fiK regime, which 
corresponds to a kinetic energy of a few kHz, the 
trapping frequencies of few Hz to kHz and potentially 
to the chemical potential of a Bosc-Einstcin condensate 
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FIG. 13: Blockade of the Rydberg excitation due to the van 
der Waals interaction. A laser with a coupling strength 
drives the transition \g) to |r) from the ground state to the 
Rydberg state. Considering two atoms, the energies of the 
pair states \g, g) and \g, r) remain almost unshifted due to 
the small polarisability of the ground state atom. The tran- 
sition \g,r) to \r,r) however is dependent on the interparti- 
cle distance r because the energy level is shifted due to van 
der Waals interaction between the two Rydberg atoms. De- 
pendent on the saturation broadening Q.r one can define a 
blockade radius below which the doubly excited state is 
suppressed. This is true as long the coupling strength Q,r of 
the excitation lasers is smaller than the linewidth T of the 
excited state. 




FIG. 14: Resonant excitation of an ensemble of two level 
atoms. When exciting an ensemble of A'^ non interacting 
atoms with a coupling strength SI the atoms undergo coherent 
Rabi-oscillations with frequency SI. In the case interactions 
and an ensemble size smaller than the blockade radius only 
one excitation is possible and the excited collective state reads 
^ = I] Iffii52,53, Ci, ...fljv). As a consequence collec- 
tive Rabi oscillations arise with a ^/N faster oscillation and 
by a factor N reduced amplitude. 



of a few kHz. All these energies are comparable with 
the spectral width of the excited Rydberg states given 
by their lifetimes. The crucial parameter for coherent 
excitation dynamics is a sufficiently strong coupling 
strength of the exciting laser fields fi^. The energy Mlji 
determines now in combination with the interaction 
energy Cqti^ all bulk properties of the driven system 
as e.g. steady state maximum Rydberg density. By 
equating these two energies Cgn^ = Kin for example 
for the 43S state and a driving field of ri/f—lMIlz one 
can estimate a minimum distance between to Rydberg 
atoms of 3.67 /im (n^ = 2 x 10^" cm^^), a distance 
severely larger then the distance between two ground 
state atoms. In other words the excitation to Rydberg 
states is strongly blockaded and only one atom out of 
rir/n sa atoms can be excited to a Rydberg state. This 
line of arguing is equivalent to the picture shown in Fig 
1131 where the power broadening 17 of the exciting laser 
field can compensate the van der Waals repulsion up to 
a minimum distance. 

Fig. [13] gives a good illustration of the blockade 
mechanism during excitation of two atoms but is only 
correct approximately, since it neglects the many body 
aspect of the excited state as described in Fig. [TH 
Within the blockade region one excitation is distributed 



over all atoms Ni, which results in a collective state 
|e) = iV'e) = 7^E^lM5i>ff2,ff3,---,e^,...gAr,), where 
we have dropped a spatial dependent phase factor 
exp(— ifcr) of the excitation light assuming a plane wave. 
The effective coupling D, to the collective excited W-state 
shows now a V^Vb enhancement due to the collective 
dipole matrix element (f/|ec?|W^) with the ground state 
1.9) = I5iiff2,53, ■■■gNb)- This W-state is actually the only 
bright state which couples to the light field, whereas the 
other Nb — 1 states are dark states as can be easily seen 
for A^6 = 2, 3 and 4: 
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(9) 



As one can easily see that in contrast to the bright W- 
state, the representation of the dark states is not unique. 
Due to the coUective enhancement of the collective exci- 
tation one has to adjust the energy argument given above 
to a collective Rabi frequency VL = y/Nbflfi. as 



(10) 



Here we introduced a coordination number which ac- 
counts for the interaction to all next and next-next neigh- 
bors for ensembles larger than the blockade radius. As- 
suming a crystalline ordered state of the superatoms with 
a hexagonal closed packing, which has the lowest en- 
ergy for repulsive van der Waals interactions, the effec- 
tive number of neighbors \s Z = 14.4. The number of 
blockaded atoms in each blockade sphere is given by 
the ratio of the density of ground state atoms divided 
by the density of Rydberg states n/rir and with this the 
density of excited Rydberg atoms and of the collective 
Rabi-frequency can be written as 



Uriv) 



( 2h 



2/5 



yzCe 
V 2h 



n(r)2/54/^ 



(11) 
(12) 



Here we have allowed for a space dependent density 
n(r) which adopts, e.g., for a Gaussian distribution for 
a thermal sample in an harmonic trap above Tc a Gaus- 
sian distribution. The Rabi frequency has been assumed 
to be constant over the whole ensemble which can easily 
be achieved by choosing sufficiently large laser beams. 
At large radii usually the density of the atomic cloud 




FIG. 15: Rydberg excitation of an inhomogeneous atom cloud 
in the superatom picture. The size rs of the blockade radius 
scales with the ground state density as n~^^^^ (see eq. 1111) . 
When the atom number per blockade sphere approaches unity 
(rir ~ n then the collective description breaks down and all 
atoms in the dilute wings of the Gaussian density distribution 
have to be treated as individual atoms. 



drops much faster than a Gaussian distribution due to a 
steeper potential for magnetically trapped clouds or the 
finite trapping region of optical dipole traps. At a cer- 
tain radius the superatom model breaks down (iVf, < 1) 
and the evolution of the driven atoms arc given by sin- 
gle atom Rabi frequency fi/j which has to be included 
in numerical simulations by setting £7 = 57^ for < 1. 
The physical situation according to the experimental ap- 
proach is depicted in Fig [15] where the inhomogeneous 
density distribution also causes a variation of the block- 
ade radius 

We start with an inhomogeneous atomic cloud with all 
atoms in the ground state and then switch on suddenly 
the coupling light field to the Rydberg state. The evo- 
lution of the excited state fraction /a, which is the total 
number of excited Rydberg atoms divided by the number 
of ground state atoms N = J n{r)dV, follows then 



n{r) 



n(r) 



17fl(r)t/2 dV, (13) 



where we have used [TT] and [T^l The integral is only 
within the volume where the blockade is effective. Fur- 
ther out one has to include the single particle dynamics 
separately. Note, that in this picture of noninteracting 
superatoms a homogeneous sample {n = const) fnit) 
would oscillate with the collective Rabi-frequency. This 
is unlikely to happen for a thermal ensemble of atoms due 
to the random spatial normal distribution of atoms lead- 
ing to a band structure of the excited states. An ordered 
array of atoms as e.g. given by a Mott-Insulator state will 
exhibit a less complex distribution of the excited states 
with a few discrete energies and coherent Rabi-oscillatios 
may be observable. The temporal evolution of fnit) for 
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a Gaussian density distribution results in a superposi- 
tion of Rabi-oscillations with varying collective frequen- 
cies and resembles a simple saturation behavior with an 
initial linear increase with a rate R and a exponentially 
approaching saturation value fn^sat as 



fR{t)=fR,sat{^ 



(14) 



This temporal evolution agrees well with the experi- 
mental observations and also with the corresponding nu- 
merical simulations. For very short times t < fl^^ the 
interaction energy lies below the saturation broadening 
and the system evolves quadratically with ^ ^j^t^- 

In the case of a Gaussian distribution one can extract 
scaling laws for R and fn^sat with the help of equations 
[TT] and [12 as a function of the peak density tlq and the 
single atom Rabi frequency ^ln 



R 



f R^sat 



n, 



4/5^2/5 



^R 



(15) 
(16) 



Keeping these results in mind, one notices that some 
care has to be taken when investigating the strongly in- 
teracting limit given by na^ ^1. It might be tempt- 
ing to associate n with the Rydberg density n^, leading 
to the relation n^r^ ~ 1. Such a reasoning, however, 
ignores the fact that the system is driven between the 
ground state and the Rydberg state. We will show in 
the subsequent theoretical analysis that the global en- 
ergy scale dominating the dynamics is Cen^ with n be- 
ing the ground state density. Treating the driven system 
as a gas of Rydberg atoms and ignoring the ground state 
atoms misses the important fact that the number of Ryd- 
berg excitations is not a conserved quantity. This is in 
strong contrast to experimental situations in which an 
atomic gas is brought into the strongly interacting limit 
by means of a Feshbach resonance; there, the number of 
atoms is indeed conserved. Nevertheless, there do exist 
several proposals, which try to convert some part of these 
extreme interaction energies into effective interaction po- 
tential for ultracold atoms by admixin g a sma ll fra ction 
of Rydberg states to the ground state jirl - ligl . Il04| . To 
get an idea of the actual length scales of the blockade 
radius Fig. [16] depicts its dependence for different values 
of the coupling field il, the ground state density and the 
principal quantum number. 



Mean-Field Theory of Strongly Interacting 
Rydberg Gases 

The dynamics on the timescale of the experiment is 
well described by the Hamiltonian in the rotating frame 
of the laser excitation with the Rabi frequency fl and the 
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FIG. 16: Fig a) shows the dependence of blockade radius for 
the 43S state according to equation [11] for a single superatom 
(Z=l) as a function of the ground state atom density Ug and 
for various coupling strenghts flft. (red: lOOHz, green: lOkHz, 
blue: IMHz, orange: lOOMHz). In b) the blockade radius is 
plotted for S-states as a function of the principal quantum 
number for a fixed coupling strength of Qr =10kHz and vari- 
ous densities Ug (red: 10^* m~^ , green: 10^^ m"^ ,blue: 10^** 
m~^, orange: 10^" m~^). The dashed line indicates in both 
figures the region, where the interparticle distance l/'^^/n is 
smaller than the blockade radius rt, which is emphasized by 
the black line in a). 



detuning A, 



pCO p(i) 

E-* ee ee 
■ 



(17) 

wher e the counter-rotating terms have been dropped 
[lOSf . It is instructive to first look at the Hamiltonian 
in the classical limit, i.e., = 0. For negative detuning 
A the ground state is a paramagnet with all spins point- 
ing down since adding a Rydberg excitation will cost en- 
ergy. However, for positive A the many-body ground 
state will have some Rydberg excitations present. For a 
fixed fraction of Rydberg excitations fu in the system. 
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the energy will be minimal when the interaction energy 
is minimized, which is the case in a crystalline arrange- 
ment. At A = there is a continuous phase transition 
from the paramagnetic to the crystalline phase [22j. In 
the strongly interacting limit it is convenient to use the 
interaction energy Eo = Cp'riP^'^ as the global energy scale 
and express the Rabi frequency and the detuning as di- 
mensionless quantities, i.e.. 



a 



and rj 



HA 



(18) 



The most straightforward method to deal with a sys- 
tem close to a continuous phase transition is mean-field 
theory. However, we cannot simply apply conventional 
Landau theory, i.e., performing a Taylor expansion of 
the energy in the order parameter, as there are no corre- 
lations in Landau theory, which is incompatible with the 
dipole blockade. Instead, we make an explicit ansatz for 
the pair-correlation function 



.921 



■IB. 



(19) 



In the blockaded region with the distance |r| being much 
smaller than the average Rydberg spacing qr the pair 
correlation function vanishes. On the other hand, at 
large distances |r| :S> an the correlation disappears and 
consequently (72(1") = 1. The transition from a strong 
suppression to the uncorrelated regime is very sharp and 
the pair correlation function can be modeled by a step 
function 



g2{r) ^ d {\r\ - an) . 



(20) 



Then, the mean field theory containing this pair corre- 
lation function is obtained by replacing the microscopic 
interaction in the Hamiltonian (|17p by the mean interac- 
tion of the neighboring atoms 

« [P^l^fR + P^fR - /l] 52(r. - r,), (21) 

which reduces the Hamiltonian to a sum of local Hamil- 
tonians iJul at each site i. In the strongly interacting 
regime with a, 77 ^ 1, the number of atoms in the block- 
aded regime is large, i.e., a^n ^ 1. This allows us to 
replace the summation over the surrounding atoms j by 
an integral over space with a homogeneous atomic den- 
sity n, i.e.. 



Cri 



-^d ^ p—d 

-nUparj , 

p — a 



(22) 



where Ad is the surface area of the c?-dimensional unit 
sphere. Then, we obtain the Hamiltonian for the i*'^ 



/« - 




FIG. 17: Excited state fraction /r of a strongly interacting 
Rydberg gas according to mean-field theory (d = 3, p = 6). 
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(23) 



In the limit of strong blockade the correlation length ^ is 
equal to an as this will minimize the interaction energy 
for a given Rydberg fraction. The correlation function 
52 (r) thus satisfies the normalization condition 



(24) 



nfnj dr[l-.g2(r)] = l, 

which provides the relation 

an. = {Vd.fBn)-^'\ 



(25) 



with Vd being the volume of the c?-dimensional unit 
sphere. The effective Hamiltonian is equivalent to a spin 
in a magnetic field h, which possesses both a transverse 
component and a longitudinal component h^, and a 
constant energy offset h^. We can diagonalize the Hamil- 
tonian using a spin rotation, i.e., i?MF = ha''^) -f /iq with 
h ~ yJK^ + hl. Here, tr^',' = cosflai*^ -|-sin0cri'^ denotes 
the Pauli matrix in the new basis while the rotation angle 
can be expressed in terms of the magnetic field compo- 
nents as tan0 = hx/h,. The ground state of the system 
is then given by {(Jz') = ^ 1, leading to the self-consistent 
equation 

a = ^UAdV!'"-' - prif^/\ (26) 

p — a 

This equation can be cast into the form 



(27) 



from which the critical exponents 5 = p/d+1 /2 and 13 = 
d/p can be determined. Solving the mean-field equation 
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for fa leads to the equation of state shown in Fig. [TTl 
with the hmiting cases 



fR - V^/" 
2d/{2p+d) 



for 
for 



: 
0. 



(28) 
(29) 



These critical exponents can also be derived using a uni- 
versal scaling function [25| . 

So far, the mean-field analysis has been centered onto 
the ground state properties of the system. However, in 
typical experiments the Rabi frequency is not varied adi- 
abatically such that the system follows its many-body 
ground state; instead, the driving lasers are switched to 
full power instantaneously, leading to a complex relax- 
ation dynamics towards a stationary state. The resulting 
critical exponents are identical to those for the ground 
state. The phenomenological superatom model can then 
be derived by calculating the energy gap h, which for 
77 = scales like h ^ ^/TTt^l with Nf, being the number 
of blockaded atoms. 



can be reduced by averaging over many different initial 
positions of the atoms according to the same distribution 
function. 

Using our numerical simulation we can also test the 
validity of our modelling of the pair correlation function 
as a step function. As can be seen from Fig. [19] the 
assumption is qualitatively correct although deviations 
can be clearly seen. There is a pronounced peak near the 
correlation length ^ and also the normalization constraint 
of Eq. (pi)) leading to ^ = is also slightly violated. 

It is now possible to investigate how the scaling of the 
saturated Rydberg fraction /sat changes when the dimen- 
sionless parameter a is varied. For this, the value of /sat 
needs to be averaged over many random configurations 
of the atoms. Then, varying a by both by changing the 
density n and the Rabi frequency fl leads to a clear data 
collapse to single line, with the numerically obtained crit- 
ical exponents being in excellent agreement with mean- 
field results 221. 



Numerical Simulations of Strongly Interacting 
Rydberg Gases 

To get a quantitative assessment of the validity of 
mean-field results it is instructive to obtain exact nu- 
merical results on the basis of the full spin Hamiltonian 
^T7\ . For a system with N atoms the dimension of the 
full Hilbert space grows exponentially like 2^; therefore, 
the exact behavior can only be calculated for a relatively 
small number of atoms. However, the strong van der 
Waals repulsion suppresses the occupation probabilities 
of many basis states, which allows us to significantly re- 
duce the Hilbert space: for each basis state we compute 
the van der Waals energy and remove the state if its van 
der Waals energy is larger than a cutoff energy Ec ■ This 
reduction leads for = 100 to approximately 10^ rele- 
vant basis states compared to the 10'^'' basis states of the 
full Hilbert space. Convergence of this method can be 
been checked by increasing Ec- 

Within mean-field theory we have already seen that the 
critical properties are identical whether we study ground 
state properties or the relaxation to a stationary state 
after a sudden increase of the Rabi frequency. As we 
are also interested in dynamical properties we perform a 
numerical integration of the time-dependent Schrodinger 
equation with the initial state being the product state of 
all spins pointing down. We place the atoms randomly 
according to a uniform distribution into a box of unit 
volume having periodic boundary conditions. Then, the 
resulting dynamics obtained using a fourth-order Runge- 
Kutta method is shown in Fig. [181 One can see the ex- 
pected relaxation towards a stationary state, although 
for A^ = 50 atoms there are considerable finite-size effects 
present in the system, especially concerning the satura- 
tion value of the Rydberg fraction //j. However, these 



Simulations of Rotary Echo Techniques 

A rotary echo setup can be used to investigate the co- 
herence time of the system subject to internal and exter- 
nal dephasing processes. Examples of external dephas- 
ing are spontaneous emission and a finite laser linewidth, 
whereas the dephasing of the single particle coherence 
due to the van der Waals interaction acts as an internal 
dephasing mechanism. In such an echo setup the Rabi 
frequency changes sign at some time Tp and the Rydberg 
fraction at the time t > Tp is measured. In absence of 
any internal or external dephasing mechanisms the dy- 
namics would be reversed, leading to /^(t = 2rp) = 0. 
Consequently, the Rydberg fraction serves as a good in- 
dicator of the coherence properties of the system. In the 
presence of dephasing the time r can be chosen such that 
/fl(r) is at a minimum. Defining the Rydberg fraction 
for the dynamics without fiipping the Rabi frequency as 



0.1 


1 1 1 


0.08 




0.06 




0.04 


+ 


0.02 


+ 








2 4 6 8 



FIG. 18: Relaxation dynamics of the Rydberg fraction fn 
(p = 6, d = 3, A = 50, a = 1/25, 77 = 0). 
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f'j^, we can define the visibility 



fkir)-fR{r) 
fnir) + Mr) ' 



(30) 



For simplicity, we consider a resonant driving, i.e., 77 = 0. 
Then, from mean-field theory we can expect v to be 
a function depending only on the intrinsic timescale 
a~^/^T. Figure [20l shows the time evolution of a typical 
echo sequence and the scaling behavior of the visibility. 
As predicted by mean-field theory there is a data col- 
lapse when the visibility is expressed as a function of the 
intrinsic timescale a~^/^T. 

The concept of the visibility can also be used to de- 
termine a (possibly time-dependent) dephasing rate 7^, 
defined by the rate equation 



-ld{T)v{T). 



(31) 



Assuming 'ydir) to be a power law in r we can directly 
express the rate up to a constant factor as 



Id = 



logv 



(32) 
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EXPERIMENTAL RESULTS ON STRONGLY 
INTERACTING RYDBERG GASES 

Strongly interacting Rydberg are best explored in 
dense ultracold clouds confined either in magnetic traps 
or in dipole traps with various methods. But also studies 
on laser cooled ensembles at much lower densities give 
valuable insight into various interaction mechanisms and 
will be discussed below. Figure [^T] shows an overview 
over our experiments on the 43 S state to date which 
have been conducted in the strongly interacting regime. 
The focus of these experiments has been lying on the 
coherent properties of the van-der-Waals interaction and 
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FIG. 19: Pair correlation function 52 (r) in the stationary state 
[p ^ 6, d = 3, N = 50, a = 1/25, ?? = 0). The length scales 
corresponding to the cutoff energy Ec and the mean Rydberg 
spacing an are shown as vertical lines. 



FIG. 20; Numerical results for a rotary echo setup, (a) Time- 
evolution of the Rydberg fraction fa with and without chang- 
ing sign of the Rabi frequency at Tp = 0.5 fi^^. (b) Scaling 
of the visibility v for different times r and different densities. 

As predicted by mean-field theory the visibility only depends 

-2/5 
on a ' r. 



how to control them. In our first experiment on strongly 
interacting Rydberg gases [31 we studied the excitation 
dynamics into Rydberg states of magnetically trapped 
clouds close to quantum degeneracy. In this work, as 
well as in all follow-up experiments, the state of choice 
was 43S of Rubidium, which yields a purely repulsive 
van der Waals potential. By this the Rydberg atoms are 
well protected against inelastic collisions as it happens 
preferential for attractive D-states (76l - [80j . Nevertheless 
some ions are also produced during the excitation into 
the 43S state and it is necessary to apply some electric 
field (^2V/cm) to extract the ions rapidly from the 
atomic cloud. As a result we could observe in [l^ a 
coherent evolution into the strongly blockaded regime 
by means of collective excitations and explained the 
scaling behavior within the phenomenological formalism 
of strongly interacting Rydberg atoms described above. 
Since the excitation rate and the steady state population 
depend on the ground state density it is possible to 
monitor the phase transition of thermal cloud to a 
Bose-Einstein condensate due to the rapid alteration of 
the density distribution ^2|. Until now the coherent 
character of the excitation has been deduced from the 
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scaling behavior of the excitation rate as a function of 
the intensity n and the couphng strength Q.^. By apply- 
ing more complex pulse sequences as e.g. in the case of 
a rotary echo it is possible to study directly coherence 
properties as depicted in Fig. [511 This has been done for 
the 43S state |9l[ and also for various D states between 
n=40 to n=45 in the vicinity of a Forster resonance [2l| . 
Also electro-magnetically induced transparency has been 
used as a coherent spectroscopy method to investigate 
the influence of interactions by means of the 43S state 
[93} . The strongly blockaded regime is accompanied by 
a nonuniform distribution of the excited Rydberg atoms 
induced by the blockade radius. A first direct observa- 
tion of these correlations has been achieved by ionizing 
the Rydberg atoms and subsequent blow up imaging of 
the resulting ions by a fluorescence detector |l06l |. Very 
recently we have used Ramsey spectroscopy to inves- 
tigate Forster resonances in the vicinity of the 44D state. 

Various aspects of the interaction mechanisms among 
Rydberg atoms, especially Forster resonances, have been 
studied already by several groups with optically trapped 
atoms [Toil - fliot and laser cooled atoms at somewhat 
smaller densities [13, El H IZll, llTll - flTTt . The sum of 
these experiments demonstrate the strength of resonant 
interactions in various approaches for different tasks, 
in which the studied effects become noticeable as line 
broadening and shifts, population transfer, and excita- 
tion blockade, which do not necessarily require through- 
out coherence. Nevertheless, the importance of resonant 
interactions lie in their tuneability, their long range char- 
acter and their anisotropy 118l | and will certainly be in- 
cluded in various experimental settings. 



CONCLUSION 

This article reviews only a small fraction of mo dern 
Rydberg physics and other review articles 0, Ill9| 
should be consulted to obtain a more complete picture. 
The research field dedicated to interacting Rydberg 
atoms is rapidly growing and splits into self-contained 
research fields. A bi g effort is made to realize scalable 
quantum computers 12Cll . 121 and quantum simulators 
[39I |40| based on individually trapped single atoms, 
whic h mi ght also be expanded to Rydberg states of 
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Actually, ensembles of strongly interacting 
Rydberg may also be useful as a source of cold ions 
since here the random potential heating is reduced by 
the blockade mechanism [23|. Although there have been 
quite some studies on the nature of Forster resonances 
there happens to be a lack of experiments actually using 
these resonances to address a specific physical problem. 
Possible fields of applications are quantum computing 
and quantum simulation. A natural research topic will 
be the study of energy transport in a quantum network 
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FIG. 21: A strongly interacting Rydberg gas for e.g. the 33S, 
43S and 53S state in Rubidium is limited to a certain do- 
main in the parameter space determined by the density n of 
the ground state atoms, the coupling strength Q, between the 
ground state atom and the Rydberg state, the strength of the 
van-der-Waals interaction Ce and the excited state lifetime. 
The radiative lifetime limit is reached when the linewidth of 
the excited state F is equal to the collective Rabi-frequency 
rjjv and the coherent excitation dynamics gets damped out. 
This restriction in the time-domain can be tightened even 
more if a strictly frozen regime is required, where the limiting 
time is now given by the motion of the atoms and by this 
the temperature of the gas. A frozen Rydberg gas is only 
present as long the atoms do not move more than a wave- 
length of the exciting light fields during excitation. At high 
densities the simple description of a van der Waals interaction 
breaks down, due to the fact that the distance between the 
atoms becomes comparable to the physical size of the Rydberg 
atoms and the Rydberg atoms start to overlap with close by 
ground state at oms as it is the case for ultra-longrange Ryd- 
berg molecules [103| . Another limitation could be inelastic 
collis ions of ground state atoms leading to three body losses 
[l07j . However, for the 43S state this is only relevant at den- 
sities above 10^^ m~'^, which are hard to prepare in the first 
place. At reduced densities the system becomes weakly in- 
teracting, whenever the blockade radius rb is smaller than 
the distance between two atoms establishing the unblockaded 
regime. For some cases there will also be a limit in the regime 
of very high Rabi-frequencies Q. due to multi-photon ioniza- 
tion processes. In the regime of strongly interacting Rydberg 
gases represented by the white area several experiments have 
been conducted with various methods as described in more 
detail in the text. The grey shaded area represents only the 
parameter space of all our experiments with Rydberg atoms 
in the 43S state ElMHal. 



20 





n = Density of ground state atoms 

Tir = Density of Rydberg atoms 

fl,. = Coupling strength of the first excitation step 

(here the 5S-5P transition in Rubidium at 780 

nm) 

Qb = Couphng strength of the second excitation 
step (here the 5P-nS,nD transition in Rubid- 
ium at 480 nm) 

flR = Effective Rabi- frequency for the two photon 
excitation 

SI = CoUective Rabi-frequency 

rii = Blockade radius 

a = Rescaled coupling strength flji 

S = Detuning with respect to the lower transition, 
here 5S-5P 

A = Detuning with respect to the two photon tran- 
sition, here 5S-nS 
rj = Rescaled detuning A 
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FIG. 22: Echo sequence depicted with Bloch-spheres. In a 
first step the ground state atoms are excited resonantly with 
a coupling strength Q indicated by the green arrow in a). 
After a certain time r the phase of the excitation field is 
fiipped by 180° resulting in an inversion of the excitation 
dynamics indicated by the blue arrow in b). In absence of 
interactions the population reaches after the a de-excitation 
time T again the ground state. This dynamics is perturbed 
if with an increasing population of the excited state the van 
der Waals interaction between the Rydberg states induces an 
effective detuning 5vdw as shown in panel c). The effective 
Rabi-frequency \/^^^'^'r'S^dw ^^'^^ larger and points also not 
anymore along the it-axes leading to a deviation of the Bloch 
vector from the uw-plane. 



[123l | as it happens also in many biological systems \124 
The non-linearity induced by the interactions do not 
only produce interesting new states of material quantum 
matter but will also show accordingly a back-action on 
the light field. This has been observed already for laser 
cooled clouds with EIT-spectroscopy [l^, 125 in terms 
of optical non-linearities. There exist various proposals 
to use the blockade effect as an optical non-linearity on 
the single photon level to generate non- classical s tates 
of light fields as e.g. Fock states [H, [Hi, [HI [l27| 
cat states [51 1 . 



and 



GLOSSARY 

Here we list the nomenclature of some important 
quantities to avoid misunderstanding due to a varying 
notations in other publications. 
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